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Abstract 

Results are presented of a two-pion interferometry (HBT) analysis in Pb+Au colli- 
sions at 40, 80, and 158 AGeV. A detailed study of the Bertsch-Pratt HBT radius 
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parameters has been performed as function of the mean pair transverse momentum 
kt and in bins of the centrality of the collision. From these results we extract model 
dependent information about the space-time evolution of the reaction. An investi- 
gation of the effective volume of the pion emitting system provides an important 
tool to study the properties of thermal pion freeze-out. 

Key words: Two-pion correlations, HBT interferometry, source size and lifetime, 
transverse expansion, beam energy dependence of radius parameters 
PACS: 25.75.-q,25.75.Gz,25.75.Ld 



1 INTRODUCTION 



Lattice QCD calculations predict a transition from confined hadronic matter 
to a state of deconfined quarks and gluons at a critical energy density e c around 
0.6 GeV/fm 3 [1]. Such energy densities are believed to be reached at CERN 
SPS energies and indeed evidence for a new state of matter created in Pb- 
induced reactions at the SPS has been announced recently [2]. It is expected 
that the reaction dynamics could differ considerably if energy densities close 
to or above e c are reached, and that the existence of a phase transition could 
manifest itself in the space-time evolution of the system [3,4,5,6,7,8,9]. In 
particular, the beam energy dependence of quantities characterizing the space- 
time evolution may give hints to the onset of deconfinement in heavy ion 
collisions. In addition, the beam energy range available at the SPS provides 
an important link between the nucleon dominated domain at the AGS and 
the pionic regime at collider energies, leading to a better understanding of the 
mechanisms relevant for thermal freeze-out. 

While single particle momentum distributions are also determined by the 
space-time evolution, the lifetime and the spatial extent of the system as well 
as the existence of collective velocity fields at the time of thermal freeze-out 
can be disentangled by the study of Bose-Einstein (BE) momentum correla- 
tions of identical pions via Hanbury-Brown and Twiss interferometry (HBT). 
The width of the correlation peak at vanishing relative momenta reflects the 
so-called length of homogeneity of the pion emitting source. Only in static 
sources can the length of homogeneity, in the following also called 'source 
radius', be interpreted as the true geometrical size of the system. In a dy- 
namic system, the occurence of space-momentum correlations of the emitted 
particles due to collective expansion generally leads to a reduction of the ob- 
served source radii. The degree of reduction depends on the gradients of the 
collective expansion velocity and the thermal velocity of the pions at ther- 
mal freeze-out. A differential analysis of the HBT correlations in bins of the 
pair transverse momentum thus provides valuable information not only on 
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the spatial extent but also on the properties of the collective expansion of the 
system [10,11,12,13]. 



2 EXPERIMENT 



The CERES/NA45 spectrometer at the CERN SPS is optimized for the mea- 
surement of low-mass electron pairs in the pseudorapidity range 2.1<i]< 2.65. 
It consists of two Silicon Drift Detectors (SDDs), located about 12 cm down- 
stream of the segmented Au-target, and two Ring Imaging Cherenkov counters 
for electron identification. To improve the mass resolution, the experimental 
setup was upgraded in 1998 by the addition of a cylindrical Time Projection 
Chamber (TPC) with radial electric drift field [14]. The TPC is located be- 
hind the existing spectrometer, 3.8 m downstream of the target. It is operated 
inside the field of a new magnet system. Analysis of the track curvature leads 
to a determination of the momentum of charged particles. The measurement 
of up to 20 space points per charged particle track allows the determination of 
the specific energy loss dE/dx with a precision of 10.5%, thereby supplement- 
ing the particle identification capability of the spectrometer. In addition to 
improving the invariant mass resolution of electron pairs, the TPC also adds 
sensitivity to hadronic observables close to midrapidity [15,16]. 



3 DATA ANALYSIS 



In this paper we present the results of a two-pion interferometry analysis in 
40, 80, and 158 AGeV Pb+Au collisions [17]. The data were taken during 
the 1999 (40 AGeV) and 2000 (80 and 158 AGeV) Pb-beam periods. In 1999, 
the segmented Au-target consisted of 8 subtargets, separated by 3.1 mm in 
beam direction, and a thickness of 25/z each, adding to a total hadronic in- 
teraction length of 0.82%. For the 2000 Pb-beam period, the target system 
was replaced by 13 subtargets of the same thickness and 2 mm spacing, with 
a total interaction length of 1.33%. An online centrality selection of approxi- 
mately the upper 30% of the total geometric cross section was applied, using 
the pulse height deposited by charged particles in the SDDs (in 1999) or in a 
Scintillator Multiplicity Counter (in 2000). The same detector information is 
used for an offline characterization of the centrality of the events measured in 
the 2000 run period, while for the 1999 data set the number of reconstructed 
charged particle tracks in the SDDs is used. Contributions from non-target 
interactions are found to be negligible. The data are presented in four bins of 
centrality, and the corresponding average numbers of participants are derived 
from a geometric nuclear overlap model [19] with o"nn=30 mb, resulting in a 
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total cross section of <tg=6.94 barn. Information about the four centrality bins 
is collected in Table 1. 



Table 1 

Definition of centrality bins. 



Centrality bin 


oV^geo 


^min ^max 


(iV part > 


1 


>15% 


5.8-6.5 


202 


2 


10-15% 


4.7-5.8 


236 


3 


5-10% 


3.3-4.7 


287 


4 


< 5% 


< 3.3 


347 



For the determination of the correlation function we combine like-sign pairs 
of charged particles detected in the TPC. A momentum dependent cut on 
the specific energy loss dE/dx in the TPC cleans the track sample from soft 
electrons and protons (see Fig. 1). Kaons are removed only for momenta below 
0.5 GeV/c. Requiring a match of the TPC track to the SDD system suppresses 
the contribution from secondary decay products and results in an increase of 
the correlation strength A but does not alter the extracted source radii. 
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Fig. 1. The specific energy loss dE/dx in the CERES-TPC in units of ADC counts 
for positive (a) and negative (b) particles as function of the particle momentum. 
The cuts used to enrich the pion sample are indicated by the dashed lines. Shown as 
solid lines are parametrizations [18] of the specific energy loss for different particle 
species according to the Bethe-Bloch formula. 



Similarly, taking the direction of the track from the SDDs rather than from 
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the TPC, which reduces the effect of multiple scattering, leaves the extracted 
source radii unchanged. Since SDD information is not available for the full data 
sample we used SDD information only for the tests just mentioned. For the 
results presented later we restricted the analysis to TPC information only. At 
the present stage of the spectrometer calibration the momentum resolution is 
5p/p=3-5% in the momentum range of interest, 0.25< p <3 GeV/c, consistent 
with the observed widths of reconstructed neutral strange hadrons A and 
K° s [17]. 

For each particle pair we calculate the mean transverse momentum 

h = ^\pt,i+M (1) 



and the pair rapidity 



1 Ei + E 2 + p z ,i + p Z)2 
= 2 ln E l + E 2 - P ^- P „ 2 



(2) 



The rapidity y m id (midrapidity) of the cm. system of the collision changes 
as function of the beam energy from y mid =2.23 at 40 AGeV to y mi d=2.91 
at 158 AGeV, while the acceptance of the spectrometer is fixed. For the 
determination of the correlation functions we restrict the analysis to pairs 
with kt > 0.05 GeV/c and -0.25< y^-y miA <0.25 (40 AGeV), -0.5< y^- 
Vmid <0 (80 AGeV), and -1< y nn -y mid <-0.5 (158 AGeV). We split the three- 
momentum difference vector q of two like-sign particles into components, 
<f = (<7iong) Qside, Qout) , to obtain detailed information about the space-time evo- 
lution. Following Bertsch and Pratt [4,5], q\ ong is the momentum difference 
along the beam direction, calculated in the longitudinal rest frame (LCMS) of 
the pair, q out is parallel to the pair transverse momentum vector k t and g S id c 
is perpendicular to q\ ong and q ont - Note that the pion mass is assumed for the 
calculation of the LCMS. 

The correlation function is defined as the ratio C 2 (q)=A 2 ( y q) / B 2 (q) where the 
'signal' A 2 (q) is the probability to find a pair with momentum difference q in 
a given event and the 'background' B 2 (q) is the corresponding mixed-event 
distribution. For the construction of B 2 (q), only tracks from events of the 
same centrality bin are combined. The number of pairs in B 2 (q) was chosen 
to be about ten times larger than in A 2 (q) to keep the statistical error of the 
background negligible. The normalization is part of the fitting procedure (see 
below). 

A separate analysis of positive and negative like-sign pairs gave within statisti- 
cal errors the same results for the extracted source radii. To improve statistics, 
we have thus combined the results of positive and negative like-sign pairs by 
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Fig. 2. Relative two track reconstruction efficiency as function of the opening an- 
gle a. For more details see text. 

adding the signal and the background distributions, respectively, before the 
construction of C 2 (q). This procedure is appropriate if the normalization con- 
stants of the positive and negative pair samples are identical, which has been 
verified. 



Due to the correlation between the momentum difference of a pair and its spa- 
tial distance in the detector, the finite two track separation power of the TPC 
can lead to a systematic bias of the correlation function. Fig. 2 shows the open- 
ing angle distribution of like-sign particle pairs divided by the corresponding 
mixed-event distribution. The resulting pair reconstruction efficiency has been 
normalized to unity at large opening angles. The pair efficiency is independent 
of the opening angle a except for very small a. We omitted pairs with open- 
ing angles a<8 mrad both in the construction of A 2 (q) and B 2 (q) to avoid 
a distortion of the correlation function due to the finite two-track separation 
power. This opening angle corresponds to a spatial distance of about 3 cm at 
the entrance of the TPC. Our studies have shown that the extracted source 
radii change only within their statistical errors if the two-track separation cut 
is varied between 6 and 10 mrad. 



Correlations in momentum space are caused not only by quantum statistics 
but also by final state interactions, the most important one being the mutual 
Coulomb repulsion. Since this is a two-body process, it depends only on the rel- 
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Fig. 3. Coulomb correlation function -F cou i as function of c/i nv before and after smear- 
ing with the momentum resolution. 

ative momentum of the pair in its rest frame qmv =j \f\(pi — P2) 2 — {E\ — E 2 ) 2 \. 
To account for this, we parametrized the qw-dependent Coulomb correlation 
function -F cou i(g inv ) derived by integration of the Coulomb wave functions over 
a Gaussian source of 5 fm size [20]. In the next step, the impact of the fi- 
nite detector resolution of <r(gi nv )=5-10 MeV/c (see below) on F cou \{q- lIW ) was 
evaluated using a Monte-Carlo simulation. The resulting Coulomb correlation 
functions -F cou i(g inv ) before and after momentum smearing are shown in Fig. 3. 

In most of the previous HBT analyses the Coulomb correction was carried out 
by applying qw-dependent weights to each pair in the background distribution 
B 2 (q). Afterwards, the Coulomb-corrected correlation function C2(q)=A 2 (q) / B 2 , C orr(q) 
was normalized in the region of large q and fit by a Gaussian parametriza- 
tion, which then should describe the pure BE correlation peak. This approach 
assumes that all particle pairs in the signal distribution are subject to the 
Coulomb correlation, which is in general not true under realistic physical and 
experimental conditions. In the following, we argue that physics and detector 
related effects lead to a reduction of the observed BE correlation strength, usu- 
ally expressed in terms of the parameter A, and that the same effects also lead 
to a similar reduction of the Coulomb repulsion. Therefore, the strength of the 
Coulomb correction applied to the data should be linked to the experimentally 
observed A parameter. 
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A reduction of the observed BE correlation from its theoretical value of two at 
vanishing relative momentum can be caused by coherent pion production, pi- 
ons from long-lived resonances and weak decays, particle misidentification, and 
finite momentum resolution. We assume that the contribution from coherent 
pion production is negligible because the observed pion phase space densities 
at SPS energies [21] are too low to expect a significant degree of coherence 
of pion production. Pions from long-lived resonances such as the a; (782) add 
a tail to the density distribution of the pion source, which then results in a 
narrow spike of the correlation function at very small relative momentum. 
Due to the finite momentum resolution of the detector the spike cannot be 
resolved. This leads to a reduction of the observed A parameter. Since the 
effective source size of pions emerging from these resonances is large, their 
mutual Coulomb interaction can be assumed to be negligible. The same argu- 
ments apply for pions from weak decays of e.g. A and K°. A contamination 
of the track sample by particles other than pions leads to pairs of like-sign 
non-identical particles. They obviously do not contribute to the BE correla- 
tion, but they experience Coulomb repulsion. However, due to their different 
masses, the strongest Coulomb correlation for those pairs does in general not 
occur at q=0 in the LCMS, which is calculated assuming the pion mass, but at 
finite q away from the pion BE peak. We have evaluated this for the CERES 
acceptance and Kir and pn pairs. In any case, we conclude that the contribu- 
tion to Coulomb correlation of non-identical pairs is weak, since the measured 
correlation functions are flat around the corresponding qi nv . Note that pairs 
containing electrons or positrons do not contribute to Coulomb correlations 
because most of them stem from decays or conversions far away from the in- 
teraction region. Pairs of identical particles which are not pions can give rise 
to BE and Coulomb correlations. However, a detailed comparison of the corre- 
lation functions of positive and negative like-sign pairs showed no significant 
difference of the extracted source radii and of the shape of the correlation 
functions at large q, although the level of contamination by non-pion pairs is 
very different in the positive and negative pair sample (see Fig. 1). Therefore 
we conclude that the contribution of non-pion pairs to Coulomb repulsion and 
BE correlation is negligible. 

In essence, we argue that only those pairs which contribute to the measured 
BE correlation strength also contribute to Coulomb repulsion. Therefore, the 
applied Coulomb correction should be scaled according to the observed A pa- 
rameter. Since the value of the A parameter is a priori unknown, we construct 
the correlation function from the signal and background distributions without 
Coulomb correction, and include the effect of Coulomb repulsion into the fit 
function: 



C 2 {q)=A 2 {q)/B 2 {q) = N ■ 1 + A' ((1 + G) • F* — 1) , (3) 
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with 



G — exp( -Rf ong 9iong -Rside^side 

— -^out^out — 2-R out i ong q , outQ , iong)) and (4) 
F* = w(k t )-(F coul (q- my )-l) + l. (5) 

The fit function contains the usual Gaussian parametrization G, with the 
Gaussian source radii i?i on g) -Rsidc, and R ut, and the cross-term R^ ut long which 
appears as a consequence of space-time correlations in non-boost-invariant 
systems [22]. The normalization N takes into account the different statistics 
used for the construction of A 2 (q) and B 2 (q). The Coulomb term F* contains 
the two-pion Coulomb correlation function F cou i(q± nv ) which was derived as 
described before. Its application requires supplementary information about the 
average q inv in each of the 10 MeV/c (gi ong , Qside, <?out)-bins, which was stored in 
separate arrays of the same bin size. Following the above argumentation, the 
strength of the Coulomb correction needs to be adjusted to the A parameter. 
At this point, an additional complication arises due to the finite momentum 
resolution, which leads to a k t dependent reduction of A by 13% to 45%, 
i.e. A=A •w(k t ). The k t dependent correction factors w{k t ) were determined 
by a Monte Carlo simulation and range from 1.15 in the lowest k t -hm. to 
1.82 in the highest k t -hm.. The inclusion of w(k t ) into the Coulomb term 
F* accounts for the apparent depletion of A due to the finite momentum 
resolution and assures the appropriate scaling of the Coulomb correction. We 
note that Eq. (3) is identical to the parametrization proposed in [23] for perfect 
resolution, e.g. w(k t )—l. 

In Fig. 4 are shown the parameters A' and \=\' -w(k t ) for 7r + 7r + and -k~it~ 
correlations at 158 AGeV (a/a geo <30%). The A parameter for negative pions 
increases with k t as expected from a decreasing contribution from electrons and 
pions from resonance decays. In contrast, for positive pions the A parameter 
decreases with k t because of an increasing contamination by protons. This be- 
haviour is in quantitative agreement with the contribution from misidentified 
particles and pions from resonance decays as seen in Monte-Carlo simulations. 

In Fig. 5 we present examples of one-dimensional projections of the corre- 
lation function, integrated over the range |q , |<30 MeV/c in the other two 
components. Also shown are the corresponding projections of the fit function 
Eq. (3). We want to emphasize that the depletion of the correlation function 
due to Coulomb repulsion is not only visible at the smallest q but also in the 
region outside the BE correlation peak. For kinematical reasons this effect 
is most pronounced in g out and increases with k t . An incorrect treatment of 
the Coulomb repulsion strength therefore affects -R S id e and R out in a differ- 
ent way, in particular leading to a /^-dependent bias of the i? ou t/-R S ide ratio. 
This is demonstrated in Fig. 6, where the source radii and the ratio -R ou t/-Rsidc 
obtained for different Coulomb strengths from central 158 AGeV events are 
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Fig. 4. Parameters A' (open symbols) and A = w(kt) ■ A' (full symbols) as function of 
kt obtained from an analysis of tt + tt + and 7r - 7r~ pairs in central 158 AGeV Pb+Au 
collisions. For explanation see text. 




Fig. 5. Projections of the two-pion correlation function in central 158 AGeV 
Pb+Au collisions (0.15<fci<0.25 GeV/c). The data are integrated over the range 
|c/|<0.03 GeV/c in the non-projected coordinates. Also shown are the results from 
a fit using Eq. (3). 
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Fig. 6. Radius parameters (panels (a-c)) and ratio R ut/Rside (panel (d)) as function 
of kt assuming different strength of the Coulomb repulsion (see text). Shown are the 
results for (1.) full Coulomb strength (open triangles), (2.) no Coulomb repulsion 
(open squares), (3.) Coulomb repulsion coupled to the A parameter as in Eq. (3) 
(full circles), and (4.) same as (3.) but with w(kt) = 1 (open circles). 

shown. We compare the results assuming full Coulomb strength (open trian- 
gles), where all background pairs were corrected by -F cou i(g inv ), with different 
implementations of the Coulomb term F* in Eq. (3). We considered the cases 
F*=l, corresponding to no Coulomb correction (open squares), F*=F cou \(q inv ) 
(open circles) as proposed in [23], and finally F*=w(k t ) ■ (F coul (g inv )-l)+l, as 
given in Eq. (5) (full circles). The parameters -Riong and -R S id e depend very 
little on the assumed strength of the Coulomb repulsion, however, the results 
for R out and correspondingly -R ou t/-Rside are very sensitive to the procedure 
employed Q. We also note that the inclusion of w{k t ) has only little effect on 
the extracted source radii. In this context it should be mentioned that our fit 
function Eqs. (3)-(5) describes very well not only the Gaussian BE correlation 
peak but also the observed undershoot of the baseline at large g out (see Fig. 5, 
right panel). 

The finite angular and momentum resolution of the TPC cause a systematic 
reduction of the extracted source radii. Corrections were determined by a 
Monte-Carlo simulation and depend on kt as well as on the extracted source 
radii themselves. The resolution in the components of q is shown in Fig. 7. 
Since a(gi ong ) and a(g S ide) are 3-7 MeV/c only, the corrections to be applied to 
-Ri on g and -R S id e are below 2%. This is demonstrated in Fig. 8 where the source 
radii from central 158 AGeV events before and after correction for the finite 
momentum resolution are shown. However, a(q ont ) deteriorates significantly as 
function of k t , which implies corrections up to 60% at high k t and consequently 
larger systematic uncertainties in R out at these momenta. 

1 Analysis of the RHIC HBT data [24] using the procedure described above should 
lead to an increase of R ou t an d R ou t /Rside [25] . Whether this can explain part or all 
of the RHIC HBT 'puzzle' remains to be investigated. 



11 



— 0.06 
o 

> 

CD 

O 0.05 



0.04 



0.03 
0.02 
0.01 




o (q in v) 



. ^long) 
" ^side) 



0.1 0.2 0.3 



0.4 0.5 0.6 

k, (GeV/c) 



Fig. 7. Resolution in (fr nv , q\ ong , q s idei and g Q ut as determined by a Monte Carlo 
simulation. 



E 
DC 



8 
6 





(b) Rside 


(C) Rout 
















O *\ 








o 




" o before correction 






o 


• after correction 









0.25 0.5 



0.25 0.5 0.25 0.5 

kt (GeV/c) 



Fig. 8. Left to right: The radius parameters R\ ng, Rside, and R ont as function of 
kt- The open symbols are before, the full symbols are after correction for the finite 
momentum resolution. 

We conclude that the main sources of systematic errors are the remaining un- 
certainties in the treatment of the Coulomb repulsion and the corrections for 
the finite momentum resolution. Varying the Coulomb strength and the mo- 
mentum resolution within realistic limits, we estimated the systematic errors 
in -Ri on g and -R S idc to be 4%. In R out , the systematic error is larger and depends 
on k t . We estimated that the systematic error in R out is 5% in the lowest and 
20% in the highest /c r bin. In the following we present the results with only 
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Fig. 9. The longitudinal radius parameter i?i ong as function of k t in different bins of 
centrality at 40, 80, and 158 AGeV. Also shown as dashed lines are fits to the data 
assuming longitudinal boost invariance (see text). The shaded regions indicate the 
systematic error. 

statistical errors at the data points. The systematic errors are indicated by 
shaded bands in the figures. 



4 RESULTS 



We present, in Fig. 9, the results for the longitudinal source radius -Ri ong as 
function of k t and in different bins of centrality for all three beam energies. At 
all beam energies and centralities under investigation the -Riong values exhibit 
a very pronounced fct-dependence. This is characteristic for the case of a strong 
longitudinal expansion, where the length of homogeneity is entirely saturated 
by the thermal velocity scale ^Tf/m t [11,12]. For the limiting case of a boost- 
invariant expansion in longitudinal direction the following relation was given 
by Sinyukov [12], which allows to extract the duration Tj of the expansion 
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Fig. 10. Centrality and beam energy dependence of the duration Tf of expansion, 
the geometric transverse source size R gco and the average transverse flow velocity 
(vt) derived from fits to the fej-dependence of i?i on g and -Rside (see text). The dashed 
lines in the right panel indicate the results obtained at 158 AGeV for a change of 
the freeze-out temperature Tf by ±20 MeV. 

(the 'lifetime') of the system: 

Rlang=Tf(Tf/mt)*- (6) 



In real sources boost-invariance is violated by the finite extension of the system 
in longitudinal direction. Sophisticated model analyses which consider more 
realistic density profiles do, however, extract lifetimes which are quite similar 
to those obtained from Eq. (6) (see e.g. [26]). Fitting expression (6) to the 
data and assuming a constant thermal freeze-out temperature Tf =120 MeV 
results in a smooth increase of ry both with centrality and beam energy from 
about 6 to 8 fm/c as shown in the left panel of Fig. 10 and Table 2. The 
choice of Tf =120 MeV is motivated by previous analyses of AGS and SPS 
data [27,28], which lead to a good description of single particle spectra and 
HBT correlations. Since the temperature enters only to the power 1/2, unre- 
alistically large changes in Tf are required to significantly affect Tf. We note 
that the increase of -Ri on g with centrality and beam energy could alternatively 
be explained by a constant duration of expansion Tf, if Tf increases by 25% 
with centrality and by 50% from the lowest to the highest beam energy. 

The results for i? S id c are shown in Fig. 11. They exhibit a weaker but still 
significant /^-dependence, indicating collective expansion also in transverse 
direction. As expected from collision geometry, -R S ide slightly increases with 
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Fig. 11. The sideward radius parameter i? s id e as function of kt in different bins of 
centrality at 40, 80, and 158 AGeV. Also shown are fits to the data (see text). The 
shaded regions indicate the systematic error. 

centrality. Comparing different beam energies, we observe that, at 40 AGeV, 
i?side at small kt is larger and the /^-dependence is systematically steeper than 
at the higher energies, while the results at 80 and 158 AGeV are very similar 
to each other. All measured -R S ide parameters at small k t are larger by almost 
a factor of two compared to the geometrical size of the initial overlap region 
of the two colliding nuclei. 

Hydrodynamically inspired model calculations predict an interplay between 
geometrical and thermal length scales. Near midrapidity, -R S idc can be approx- 
imated by [29,30,31]: 



#side = Rgeo/ (1 + rrif r] 2 f/ T f) 2 ■ 



(7) 



This expression can be used to determine the 'true' geometric transverse size 
R geo of the source at freeze-out. We estimate from a model dependent analysis 
of the fit parameter v}/Tf the average transverse flow velocity. This fit param- 
eter contains the thermal freeze-out temperature Tf and the transverse flow 
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Fig. 12. Panels (a-b) show the la- and 30<7- contours of the fit of Eq. (7) to i? s ide(&t) 
from central 40 AGeV and 158 AGeV collisions in the Rg eo -v} /I/-plane. Also in- 
dicated by the dotted lines is the location of x 2 -valleys, and the position of the 
minimum by the crosses. In panel (c) are shown the x 2 - varues computed along the 
valleys indicated in (a-b) and plotted as function of rq/Tf. 



rapidity rjf, which characterizes the strength of transverse flow for a given flow 
profile and transverse density distribution. We applied a fit of this expression 
to the data, as indicated by the dashed lines in Fig. 11. The fits describe the 
data very well at all centralities and energies. However, we achieve reason- 
able results and errors for the fit parameters only at 80 and 158 AGeV, as 
shown in Table 2. A closer inspection of the x 2 -contours of the fits, shown in 
Fig. 12, shows that the parameters R geo and rjj/Tf are highly correlated and 
that the x 2_con tours are very flat and without discernible minimum for the 
40 AGeV data. However, the \ 2 values indicate that the 40 AGeV data are 
also consistent with fit parameters similar to those obtained for the higher 
energies, at in fact the same level of confidence as at the higher energies. We 
conclude that a separation between source geometry and flow dynamics using 
expression Eq. (7) is not possible because of the steepness of R S id c (k t ) at 40 
AGeV. 



The results obtained at 80 and 158 AGeV are very similar to each other: We 
observe a smooth increase of R gco with centrality from about 5.5 to 7 fm, as 
shown in the middle panel of Fig. 10. Assuming a thermal freeze-out tempera- 
ture of T/=120 MeV and a box-shaped density distribution, the average trans- 
verse flow velocity (v t ) can be derived from the fit parameter i]}/Tf [32,26], 
as shown in the right panel of Fig. 10. For the most central events at 80 and 
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158 AGeV we obtain (vt) ~ 0.55c and very little centrality dependence. A 
variation of the assumed freeze-out temperature Tf by ±20 MeV alters the 
extracted flow velocities only within their statistical errors, as demonstrated 
for the 158 AGeV results by the dashed band in Fig. 10 (right panel). These 
findings are similar to results from previous analyses [28,26,33]. 



Table 2 

Summary of freeze-out parameters obtained from fits to the ^-dependencies of i?i on , 
and -Rgidc (see text). 





Cent, bin 


Tf (fm/c) 


.Rgeo (fm) 


rij/Ti (GeV- 1 ) 


(vt)/c 


40 AGeV 


1 


5.95±0.14 


6.oii^ 5 7 


3-74±f:|i 


- 


2 


6.34±0.13 


10.101"? 


15.43l 3 7 5 4 f 


- 


3 


6.49±0.11 




12-84li 1 9 f 




4 


6.63±0.10 


12.08l™j 


19.06+^ 8 7 f 




80 AGeV 


1 


6.49±0.18 


k qi+1.04 


o crq+2.60 


0-461°:^ 


2 


6.98±0.10 


r- no+0-57 
d - yz -0.43 


q rv^i+1.33 


o ™+ - 07 

u - ou -0.07 


3 


7.29±0.09 


6 23+ 52 


o nq+1-14 


^n +0 - 06 

U - OU -0.06 


4 


7.30±0.09 


7 iq+0.78 


4.5iii;l 


n c;o+0.07 
U - b8 -0.06 


158 AGeV 


1 


7.34±0.09 


r cq+0.44 


2 9fi +L11 
z - yD -0.82 


n 4q +0.06 


2 


7.46±0.08 


r qn+0.44 
°- yu -0.36 


q cq+1.13 


r, c;q+0.06 
u - oo -0.05 


3 


7.87±0.06 


5-8618:* 


9 zK+0.56 
z - 4O -0.47 


46 +0 - 04 


4 


8.23±0.05 


6-8618:1? 


q qi+0.70 


u - oo -0.03 



It has been proposed that the existence of a strong first order phase transistion 
and an accordingly long lived mixed phase would be observable by a large 
outward radius R out compared to -R S id e , indicating a long duration of pion 
emission At [4,5,6,7,8,9]: 



Ar2 = ~M ( R lnt ~ R lide), ( 8 ) 
Pt 

where f3t~kt/mt is the average transverse velocity of the pion pair in the LCMS 
(/3/~0). The dependence of the R out parameter on k t is shown in Fig. 13. For 
comparison, -R S id c is also shown by the solid lines. We observe a smooth increase 
of R out with centrality and beam energy. At 158 AGeV, the ratio -R ou t/-Rside 
is larger than unity at small k t . While -R out /-R S idc<l was observed at large 
k t recently at RHIC, we see at the highest beam energy a ratio consistent 
with one, within the systematic uncertainties which are mainly imposed by 
the large corrections to -R ou t- 



17 




Fig. 13. The outward radius parameter R ont as function of k t in different bins of 
centrality at 40, 80, and 158 AGeV. The solid line represents the corresponding 
results for i? s id e - The shaded regions indicate the systematic error. 

The data do not support the scenario of a long-lived source at SPS: All ob- 
served radii R out are similar to -R S ide at all beam energies and centralities. At 
the lowest beam energy i? S id e is even slightly larger than i? ou t- We note that 
i? out <-R S idc w as suggested as an indication for sources with surface dominated 
emission [34], such as emission from an expanding shell. At 158 AGeV, the 
data are consistent with a short but finite emission duration of about 2-3 fm/c, 
in agreement with results reported by NA49 previously [28]. 

The existence of a non- vanishing cross term R^ ut long would indicate space-time 
correlations of the pion emission points [22]. Such correlations occur if the 
longitudinal flow velocity of the source element with highest emissivity does 
not coincide with the average longitudinal velocity of the pions. For symmetry 
reasons, i?out,iong * s expected to vanish at midrapidity and in sources with 
longitudinal boost-invariant expansion. Our results for the cross term R^ ut long 
are shown in Fig. 14. While R^ ut long is consistent with zero at 40 AGeV, we 
observe a smooth increase of its absolute value with beam energy at small k t . 
The small but non-zero results at the higher beam energies may be caused 
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Fig. 14. The cross terms R^ ut long as function of k t in different bins of centrality at 
40, 80, and 158 AGeV. 

by a slight shift of the rapidity acceptance of the spectrometer away from 
midrapidity as the beam energy increases. 

Pion interferometry data published by experiments at AGS and RHIC together 
with the data presented here can be combined to perform a systematic study 
of the source parameters over a wide range of beam energies. In Fig. 15 are 
shown the /cj-dependences of R\ ong , R S ide, and R out in central Pb(Au)+Pb(Au) 
collisions near midrapidity [35,24]. No dramatic variation of the source param- 
eters can be observed. However, a closer inspection reveals interesting features. 
The parameter -Ri on g is approximately constant from AGS to the lower SPS 
energies, but starts to increase significantly within the SPS regime and to- 
wards RHIC, indicating a smooth increase of the lifetime. -R S ide is gradually 
decreasing at small k t up to top SPS energy, connected with a continuous 
flattening of the /^-dependence. At RHIC, -R S id e is again larger than at the 
SPS while the shape is not yet well measured. In fact, the -k~tt~ sample in- 
dicates further flattening. Naively, the flattening would indicate a decrease of 
the radial flow velocity as function of beam energy. This is in contradiction 
to the present interpretation of single particle transverse mass spectra, which 
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Fig. 15. Compilation of HBT radius parameters near midrapidity in central 
Pb(Au)+Pb(Au) collisions at AGS [35], SPS (present data), and RHIC [24] en- 
ergies. 

indicate an increase of radial flow with beam energy [36]. The parameter R ont 
shows a rather weak energy dependence and a slight minimum around the 
lowest SPS energy, where we find i? out /i? side <l. 

An investigation of the freeze-out conditions can be performed by relating the 
measured source parameters to an effective freeze-out volume: 

V f = (27r)§it> longJ R s 2 ide , (9) 

computed for Gaussian density distributions in all three spatial dimensions. 
In the presence of collective expansion, Vf does not comprise the total volume 
of the pion source at freeze-out but rather reflects a volume of homogeneity. If 
thermal freeze-out would occur at constant density [37], a linear scaling of Vf 
with the charged particle multiplicity would be expected. Fig. 16 shows V), 
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Fig. 16. The freeze-out volume Vf calculated at 0.15<fct<0.25 GeV/c (see text) as 
function of the centrality expressed here in terms of the number of participants at 
40, 80, and 158 AGeV. 

computed at 0.15<A;i<0.25 GeV/c, as function of the number of participants 
at 40, 80, and 158 AGeV. An approximately linear scaling with N paxt is in- 
deed observed at all three energies. Since the number of charged particles was 
found to scale close to linear with N paxt at SPS [38,39,40], this is consistent 
with the assumption of a constant freeze-out density at a given beam energy. 
However, the observed beam energy dependence is at first glance surprising: 
the increase of pion multiplicity at midrapidity by about 70% from 40 AGeV 
to 158 AGeV [41] is not reflected in a corresponding increase of Vf. Obvi- 
ously, the freeze-out volume scales with multiplicity as long as multiplicity 
is controlled via centrality, but it does not scale accordingly as multiplicity 
changes with beam energy. We note that, according to experimental results, 
the ratios of particle species change very little with centrality for not too pe- 
ripheral collisions, whereas they change significantly with beam energy. In a 
separate publication [42] we show that thermal pion freeze-out in general does 
not occur at constant density pf—N/Vf, if relative abundancies of hadrons 
change. This is due to the different cross sections of pions with e.g. pions and 
protons, which explains that freeze-out occurs rather at constant mean free 
path \f=(pf(j e f[)~ 1 , where cr cS is the effective cross-section of the system with 
pions, averaged over all particle species in the final state: a e s=l/N ■ J2i Ni<j ni . 
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5 SUMMARY AND CONCLUSIONS 



We have presented a systematic study of two-pion interferometry data at SPS 
energies. We observe a smooth centrality and beam energy dependence of all 
HBT parameters; no indications for drastic changes of the space-time evolu- 
tion have been found. The data at all centralities and beam energies under 
investigation suggest the existence of a pion source which rapidly expands in 
longitudinal and transverse direction. Assuming longitudinal boost-invariance, 
a duration of expansion of about 6-8 fm/c can be determined from the pair 
transverse momentum dependence of the longitudinal radius parameter i?i on g, 
smoothly increasing with centrality and beam energy. The expansion is fol- 
lowed by a rather sudden freeze-out with a short duration of pion emission. At 
the time of freeze-out, the transverse r.m.s. radius of the system has increased 
by about a factor two compared to the dimension of the initial overlap region 
of the two colliding nuclei. This implies that the transverse area of the reac- 
tion has quadrupled. The transverse momentum dependence of the parameter 
.Rside indicates a high degree of collectivity and a radial flow velocity of about 
(v t )~0.55c in the most central collisions. 

Studying the beam energy dependence of i?i ong , there is a monotonic increase 
from top AGS energies to RHIC. Over the same range of beam energies there is 
a subtle but systematic flattening of the /^-dependence of -R S id e - More studies 
are needed to understand the origin of this effect. The centrality dependence 
of the effective freeze-out volume Vf is consistent with the assumption of pion 
freeze-out at constant particle density. However, this simple picture breaks 
down if different beam energies are compared. The key to understand this 
effect is in the consideration of the cross sections of different particle species 
with pions, as the relative abundances change with beam energy. 
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